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solving  the  Grad-Shaf ranov  equation  with  the  footpoint  displacement  and  the 
entropy  prescribed.  It  is  shown  that  no  multiple  solutions,  or 

bifurcations,  result.  This  is  to  be  contrasted  with  the  approach  of 
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19.  ABSTRACTS  (Continued) 

displacement  or  entropy,  as  opposed  to  or  p,  must  be  specified  are 
discussed.  It  is  argued  that  these  conditions  are  more  likely  to  occur  in 
the  corona  than  those  conditions  under  which  and  p  may  be  prescribed. 
The  lack,  of  bifurcations  indicates  that  equilibrium  will  not  be  lost  as  the 
footpoint  displacement  or  entropy  is  increased.  The  limiting  configurations 
for  infinite  footpoint  displacement  and  infinite  entropy  are  also  discussed. 
It  is  shown  that  although  the  current  density  does  become  somewhat  peaked, 
the  total  current  in  the  peak  region  decreases  as  the  system  is  sheared  (or 
heated).  In  fact,  the  current  in  this  peak  region  contains  a  rapidly 
decreasing  fraction  of  the  total  current  so  that  the  limiting  configuration 
is  not  one  in  which  the  current  is  concentrated  into  a  current  sheet. 
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EQUILIBRIUM  OF  SOLAR  CORONAL  ARCADES 


I.  Introduction 

Numerous  fundamental  issues  remain  unresolved  concerning  important 
energetic  phenomena  in  the  solar  corona.  For  example,  the  possible 
mechanisms  of  solar  flares,  coronal  mass  ejections  and  various  eruptive 
processes  are  not  well  understood.  In  these  energetic  events,  energy 
release  is  manifested  in  the  form  of  heating,  radiation  and/or  mass  motion 
in  the  coronal  plasma.  In  addition,  the  corona  evidently  undergoes 
continual  and  pervasive  (but  perhaps  less  explosive)  heating.  Various 
theories  have  been  proposed  to  explain  these  observed  phenomena  (for  an 
overview,  see,  for  example,  Sturrock  1980).  A  substantial  number  of  them 
assume  that  magnetic  fields  are  the  ultimate  source  of  energy.  This 
underlying  assumption  is  strengthened  by  the  recent  observations  (e.g.,  the 
Skylab  mission,  1973)  showing  that  the  corona  is  a  highly  complicated  system 
with  magnetic  fields  and  currents  presumed  to  control  the  structuring  and 
dynamics  of  coronal  plasmas.  As  a  result,  it  is  important  to  understand  the 
equilibrium  and  dynamical  properties  of  magnetic  structures  in  the  corona. 

In  studying  coronal  magnetic  fields,  two  dimensional  magnetohydrodynamic 
(MHD)  equilibrium  models  have  received  considerable  theoretical  attention. 
Such  models  are  important  both  because  of  their  relative  simplicity  and 
because  many  fundamental  issues  that  arise  in  2D  may  also  exist  in  three 
dimensional  systems.  For  example,  linear  arcades  may  be  relevant  to 
magnetic  structures  associated  with  neutral  lines,  where  the  normal 
component  of  the  magnetic  field  R  is  zero.  Such  structures  have  been 
inferred  from  observations.  An  arcade  may  have  a  considerable  amount  of 
magnetic  shear  but  typically  varies  weakly  in  the  direction  along  the 
neutral  line.  As  another  example,  cylindrically  symmetric  geometry  has  been 
used  to  model  sunspot  fields  (Barnes  and  Sturrock  1972;  Yang,  Sturrock,  and 
Antiochos  1986). 

One  frequently  invoked  scenario  in  which  the  magnetic  energy  of  a 
structure  in  the  corona  is  thought  to  be  built  up  and  stored  is  the 
following;  as  coronal  magnetic  field  lines  anchored  in  the  photosphere  are 
sheared,  the  structure  evolves  through  a  sequence  of  quasi-equilibrium 
configurations,  stressing  the  field  lines  and  building  up  magnetic  energy. 
As  the  twist  of  field  lines  in  the  corona  is  increased,  equilibrium  limits 
may  be  exceeded,  resulting  in  catastrophic  loss  of  equilibrium  and  sudden 
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release  of  magnetic  energy.  This  scenario  is  an  appealing  one  because  of 
the  ubiquity  of  coronal  magnetic  fields.  A  similar  scenario  in  which  loss 
of  equilibrium  occurs  as  a  result  of  increasing  the  pressure,  a  measure  of 
the  internal  energy  of  the  plasma,  has  also  been  suggested.  The  issues 
addressed  in  this  paper  pertain  to  sequences  of  equilibria  and  loss  of 
equilibrium.  We  believe  that  these  issues  are  relevant  in  both  two-  and 
three-dimensional  structures. 

In  this  paper  we  specialize  to  straight  arcade  equilibria  such  that  3/3z 
=  0  for  all  variables.  Figure  1  shows  a  schematic  drawing  of  a  linear 
arcade  field  line  configuration.  Such  arcades  have  been  studied  in  the 
context  of  quasistatic  evolution  of  equilibria  (Barnes  and  Sturrock.  1972; 
Low  1977,1982;  Jockers  1978;  Birn,  Goldstein,  and  Schindler  1978;  Heyvaerts 
et  al.  1982;  Yang,  Sturrock,  and  Antiochos  1986;  Zwingmann  1987;  Priest 
1988;  Klimchuk,  Sturrock,  and  Yang  1988)  as  well  as  dynamical  evolution  (Vu, 
Bao,  and  Tandberg-Hanssen  1987;  Mikic,  Barnes,  and  Schnack  1988.)  Figure  1 
shows  schematically  the  coordinate  system  and  a  representative  field  line. 
The  field  lines  are  labelled  by  the  flux  =  A2(x,y).  The  footpoints  are 
displaced  in  z  by  an  amount  d(<j/)  in  the  photosphere  S.  We  adopt  the  usual 
ideal  MHD  line-tying  condition  in  the  photosphere.  Then,  the  footpoint 
displac€:ment  can  be  specified  by  plasma  motion  in  the  photosphere.  Arcade 
equilibrium  studies  have  been  previously  performed  with  the  idea  that 
magnetic  energy  W  =  /dVB^/Srt  increases  as  the  shear  is  increased  and  that 
because  of  this  increase,  MHD  instability  or  loss  of  equilibrium  is 
possible.  A  standard  method  for  obtaining  equilibrium  solutions  is  to  solve 
the  Grad-Shaf ranov  equation.  The  most  natural  approach  to  solving  this 
equation  for  a  linear  arcade  is  to  specify  the  axial  component  of  the 
magnetic  field  and/or  the  pressure  p('^).  [In  the  following,  when  we 
write  p  =  PC'!*),  we  assume  that  the  relevant  length  scales  in  the  corona  are 
smaller  than  the  gravitational  scale  height  in  the  corona.]  Thus,  in  many 
of  the  earlier  papers,  loss  of  equilibrium,  or  bifurcation,  has  been  studied 
in  the  context  of  arcade  equilibria  with  B^  specified  (for  an  overview,  see, 
for  example,  Birn  and  Schindler  1981).  Another  approach  is  to  treat  the 
equilibrium  problem  as  a  time  dependent  one  involving  Clebsch  variables  and 
introduce  artificial  viscosity  or  drag  (Chodura  and  Schluter  1981;  Yang, 
Sturrock  and  Antiochos  1986;  Klimchuck,  Sturrock  and  Yang  1988).  This 
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"magneto-frictional"  method  allows  the  system  to  find  a  stable  equilibrium 
for  specified  footpoint  displacement. 

By  solving  the  Grad-Shaf ranov  equation  with  prescribed  axial  field 
bifurcations  have  been  found  in  the  two-dimensional  force-free  limit 
(Jockers  1978).  Having  found  multiple  solutions,  Jockers  then  pointed  out 
that  it  may  be  more  physically  significant  to  specify  the  footpoint 
displacement  rather  than  B^  and  that  bifurcations  may  not  exist  with  respect 
to  specification  of  footpoint  displacement  d(\J;).  The  reason  is  that,  in  the 
idealized  limit  of  perfectly  conducting  photosphere  and  corona,  the 
specified  footpoint  displacement  can  be  conserved  as  the  system  adjusts  to 
find  its  equilibrium,  while  the  coronal  plasma  posseses  no  physical 
mechanism  which  can  hold  (or  equivalently  the  transverse  current 

perpendicular  to  z)  fixed  during  this  process.  It  also  depends  upon  the 
assumption  that  inertia  dominates  magnetic  forces  in  the  photosphere  (py^/2 
>>  B^/8n)  and  therefore  the  actual  fields  obtained  in  the  arcade  equilibrium 
do  not  influence  the  footpoint  motion.  In  a  recent  comprehensive  and 
detailed  study,  Zwingmann  (1987)  obtained  sequences  of  arcade  equilibria  by 
varying  footpoint  displacement  and  pressure.  The  conclusion  of  this  work  is 
that  no  multiple  solutions  (i.e.,  bifurcations)  exist  if  footpoint 

displacement  is  prescribed,  confirming  the  suggestion  of  Jockers  (1978),  but 
that  multiple  solutions  do  exist  if  pressure  p(\J()  is  prescribed.  Thus,  it 
was  argued  that  the  onset  conditions  for  solar  eruptive  processes  might  be 
determined  by  bifucations  with  respect  to  increasing  the  pressure.  In  this 
work,  Zwingmann  used  an  iteration  method  developed  by  Keller  (1977).  Priest 
(1988)  has  also  suggested  that  if  the  pressure  is  increased  beyond  a  certain 
critical  point,  eruptive  motion  (e.g.,  coronal  mass  ejections)  may  result 
due  to  loss  of  equilibrium.  Recently,  Klimchuck,  Sturrock  and  Yang  (1988) 
also  found  that  the  magnetic  energy  of  a  linear  dipole  field  increased 
indefinitely  with  increasing  footpoint  shear. 

Clearly,  it  is  important  to  distinguish  between  specifying  B2('J')  and 
specifying  footpoint  displacement  d(<|»)  because  the  bifurcation  properties 
are  completely  different.  An  analogous  distinction  exists  between 

specifying  the  pressure  p(vt;)  and  specifying  the  entropy  s(\|;),  both 
quantities  being  measures  of  internal  energy  of  the  plasma.  The  entropy  can 
be  defined  in  terms  of  the  heat  contained  between  flux  surfaces  (Sec.  II). 
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In  the  limit  of  adiabatic  coronal  plasma,  the  entropy  is  a  conserved 
quantity.  In  the  corona,  we  expect  perpendicular  thermal  conduction  to  be 
smaller  than  parallel  conductivity  and  different  flux  surfaces  are 
effectively  thermally  isolated.  If,  in  addition,  the  parallel  heat  flux 
between  the  corona  and  the  photosphere  is  small  enough,  then  the  corona  is 
thermally  isolated  along  the  field  lines  from  the  photosphere.  If  these 
conditions  are  satisfied,  then  the  adiabatic  limit  may  be  approximately 
valid  and  the  entropy  sfij/),  like  the  footpoint  displacement,  is  conserved  as 
the  system  adjusts  to  find  its  equilibrium.  Here,  it  is  implicitly  assumed 
that  the  Alfven  speed  or  magnetosonic  speed  is  sufficiently  fast  that 
radiation  and  thermal  conduction  to  the  photosphere  are  negligible  on  the 
time  scale  of  relaxation  to  equilibrium.  We  also  discuss  in  Sec.  II  an 
alternative  limit  of  an  isothermal  model,  corresponding  to  clamping  of 
temperature  by  radiative  processes  or  by  parallel  thermal  conduction.  This 
limit  has  essentially  the  same  mathematical  form  as  the  adiabatic  limit,  but 
with  adiabatic  index  of  unity.  In  contrast,  it  appears  that  there  is  no 
such  limiting  case,  arising  from  a  reasonable  energy  equation  for  the 
coronal  plasma,  in  which  the  pressure  p('l')  is  conserved  as  the  flux  surfaces 
adjust  themselves  to  find  an  equilibrium. 

In  the  present  paper,  we  present  a  method  to  solve  the  Grad-Shafranov 
equation  for  linear  arcades  with  prescribed  footpoint  displacement  d(\|/)  or 
entropy  s(<//).  Using  this  technique,  we  consider  the  bifurcation  properties 
of  equilibrium  arcades.  We  find  that  no  multiple  solutions  exist  if  df^J/)  is 
increased,  in  agreement  with  Zwingmann  (1987).  However,  we  find  that  there 
are  no  multiple  solutions  if  the  entropy  rather  than  the  pressure,  is 
specified.  In  situations  where  the  entropy  rather  than  the  pressure  is  to 
be  specified,  the  bifurcations  with  respect  to  the  pressure  are  of 
mathematical  rather  than  physical  significance. 

We  also  consider  questions  regarding  the  formation  of  sheet  current  in 
"open  field  configurations"  which  have  been  postulated  as  the  limit  of 
infinite  footpoint  displacement  (Barnes  and  Sturrock  1972;  Yang,  Sturrock, 
and  Antiochos  1986;  Aly  1984,  1985).  This  configuration  is  a  potential 
field  (i  =  0)  except  along  the  current  sheet.  The  development  of  sheet 
currents  ("tangential  discontinuities"  of  the  magnetic  field)  has  been  an 
active  area  of  investigation  since  it  was  suggested  by  Parker  (1972,  1983) 
as  a  possible  source  of  coronal  heating.  See  also  van  Ballegooijen  (1985), 
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Zweibel  and  Li  (1987),  Low  and  Wolfson  (1988)  and  Antiochos  (1988).  We  find 
that  as  d((ti)  is  increased,  the  current  density  does  become  more  peaked 
but  that  the  maximum  value  of  increases  only  slowly,  if  at  all,  with 

footpoint  displacement,  and  that,  relative  to  the  total  current  a 
decreasing  fraction  of  current  is  carried  by  the  current  peak.  We  conclude 
that  the  field  does  not  approach  the  postulated  open  field  configuration. 
By  the  same  token,  if  the  entropy  s(i|»)  is  increased,  the  current  density 
becomes  increasingly  peaked  but  a  decreasing  fraction  of  current  is  carried 
by  the  peak  relative  to  the  total  current  I^. 

In  Sec.  II  we  formulate  the  MHD  equilibrium  problem  for  a  straight 
arcade  with  pressure  and  gravity.  We  discuss  the  specification  of  footpoint 
displacement  and  entropy  and  relate  the  former  to  the  familiar  Clebsch 
representation  for  fields.  In  Sec.  Ill,  we  discuss  the  computational  method 
used  to  obtain  the  equilibria  when  footpoint  displacement  or  entropy  is 
prescribed.  We  then  show  results  obtained  by  this  method  for  various 
profiles  of  photospheric  flux,  footpoint  displacement,  and  entropy.  We 
compare  these  results  with  some  simple  analytic  cases  and  show  estimates 
based  upon  scaling  for  large  footpoint  displacement  or  entropy.  The 
relevance  to  the  formation  of  a  sheet  current  as  d  ■*  ®  or  s  ”  is 
discussed.  Sec.  IV  contains  discussions  of  the  results  and  their 
significance  to  coronal  observations.  The  salient  results  are  summarized  in 
Sec.  V. 
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II.  Linear  Arcade  Equilibria 

a)  Basic  Equations 

The  geometry  considered  in  this  paper  is  that  of  a  linear  arcade  with  a 
straight  neutral  line.  We  choose  coordinates  as  shown  in  Fig.  1,  with  Cy 
the  unit  vector  normal  to  the  surface  S(y  =  0)  corresponding  to  the 
photosphere  and  x  labelling  the  non-ignorable  direction  in  S.  We  adopt  the 
usual  assumption  that  the  corona  has  infinite  electrical  conductivity  (on 
the  time  scales  of  interest)  and  that  the  photosphere  is  infinitely 
conducting  and  massive.  Then,  the  field  line  footpoints  are  line-tied  in 
the  surface  S.  In  such  a  geometry  where  z  is  the  ignorable  direction,  the 
magnetic  field  can  be  represented  by 


B  =  S'!*  X  §2  +  8262  ,  (la) 

where  the  flux  function  'J/(x,y)  is  the  z-component  of  the  vector  potential. 
The  current  density  i  =  2  x  B  (in  normalized  units  c/4n  =  1)  is  then  given 
by 


i  =  yBg  X  e2  -  • 


(lb) 


Including  the  effects  of  pressure  p  and  gravity,  the  equation  of  motion  is 


i  X  B  -  2p  -  mgney  =  0, 


(2) 


where  n  is  the  ion  (or  electron)  density  and  m  is  the  ion  mass.  The 
component  of  eq.  (2)  in  the  z-direction  gives  B‘2B2  =  0  which  implies 
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Bz  = 


(3) 


The  quantity  related  in  a  simple  manner  to  the  transverse 
current  per  unit  length  in  the  z  direction.  Consider  a  closed  path  around  a 
rectangle  defined  by  <  ^2»  0  <  z  <  L  on  the  photospheric  surface  (S). 
Ampere's  lav  <jB'dl  =  I  then  yields  -  B2('|'2)J  as  the  transverse 
current  per  unit  length  in  z  flowing  between  the  flux  surfaces  ij;  =  and  ♦ 
=  y^2'  Writing  the  pressure,  temperature  and  density  as  functions  of  the 
nonor thogonal  coordinates  and  y,  we  find  that  the  component  of  eq.  (2) 
parallel  to  B  leads  to 


3y 


P('J'»y) 


■"gp(*i>.y) 

2kT(<J»,y) 


(4a) 


where  k  is  Boltzmann  constant,  T  is  the  electron  temperature  and  m  is  the 
ion  mass.  This  equation  can  be  solved  to  give 


p(*,y)  -  P(*,0)exp(-  )  ■ 

For  the  special  case  T  =  T(\J/),  we  find 


(4b) 


n(*,y)  =  no(<|»)expt-  mgy/2kT(g>)  J .  (4c) 

Here,  p(^,0),  n(<J»,0)  are  the  pressure  and  density  on  S.  Equation  (4a) 
simply  states  that  the  plasma  is  in  hydrostatic  equilibrium  along  the  field 
lines.  Substituting  eqs.  (3)  and  (4b)  in  the  2'J'  component  of  eq.  (2)  leads 
to  the  Grad-Shafrarov  (G-S)  equation 


.  (5a) 

where 

dB 

•  W  ’ 

and 
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p(vj<,y)  =  2n(«i/,y)kT(i/;,y). 


(5b) 


In  the  remainder  of  this  work,  we  adopt  the  assumption  that  y  is  much 
smaller  than  the  gravitational  scale  height  in  the  corona  h  =  2kT/mg  so 
that  force  balance  along  the  field  lines  is  achieved  by  having  pressure 
constant  along  the  field  lines  in  the  corona,  giving  p  =  p(<J/).  We  will 
also  assume  that  paralle  thermal  conduction  dominates  in  the  corona  (but 
not  in  the  transition  zone,  chromosphere,  or  the  photosphere)  so  that  T  = 
T(iJ/).  It  follows  from  eq.  (5b)  that  density  is  also  constant  along  field 
lines,  n=  n^Cij;).  An  important  special  case  occurs  if  the  plasma  beta  P  = 
2p/B^  is  small,  giving  the  force  free  condition  p  =  0. 

The  boundary  conditions  which  we  specify  on  eq.  (5)  are  that  the  flux 
\l>(x,0)  is  given  on  the  surface  S  and  is  assumed  to  be  constant  (zero, 
without  loss  of  generality)  for  (x,y)  ^  This,  of  course,  determines  the 
component  By  =  -8<j/(x,0)/3x  normal  to  S.  In  practice,  these  conditions  are 
replaced  by  those  specifying  (J/(x,0)  for  a  finite  interval  0  <  x  <  a  on  S 
and  by  requiring  \l/  to  be  zero  on  the  surfaces  x  =  0,  x  =  a  and  y  =  L.  See 
Figure  2.  Ve  will  return  to  the  implications  of  these  boundary  conditions 
on  a  box  of  finite  size  (in  x,y)  in  Sec.  III.  Mathematically,  the  G-S 
equation  (5)  is  usually  posed  as  a  nonlinear  boundary  value  problem  by 
specifying  the  nonlinear  functions  and  This  specification  can 
lead  to  multiple  solutions  (bifurcations).  However,  it  was  suggested  by 
dockers  (1978)  that  B^  is  not  a  physically  specifiable  quantity  in  the 
corona  and  that  bifurcations  may  not  exist  if  the  footpoint  displacement  is 
specified.  This  point  was  recently  demonstrated  by  Zwingmann  (1987),  who 
also  found  bifurcations  with  respect  to  specifying  the  pressure  p(i</). 

b)  Specif iability  of  Physical  Quantities 

In  this  section,  we  discuss  the  specif iabili ty  of  a  number  of  physical 
quantities.  For  this  purpose,  it  is  helpful  to  use  an  imaginary  two-step 
process.  Consider  a  system  in  equilibrium  (e.g.,  i  x  B  =  0).  Suppose  a 
physical  quantity  (e.g.,  footpoint  displacement)  is  changed  holding  all 
other  quantities  fixed  (first  step).  Then  the  system  is  no  longer  in 
equilibrium  and  it  must  relax  to  a  new  equilibrium  (second  step).  If  the 
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physical  quantity  in  question  is  a  conserved  quantity  during  the  relaxation 
process,  then  it  retains  the  given  value  after  reaching  the  new 
equilibrium.  Such  a  conserved  quantity  can  be  specified.  If,  on  the  other 
hand,  the  quantity  is  not  a  conserved  quantity,  then  it  also  changes  during 
the  relaxation  process  so  that  it  cannot  be  held  fixed.  In  this  case,  a 
specified  value  need  not  be  attainable. 

For  example,  if  the  footpoint  displacement  is  increased  from  that  of  an 
initial  equilibrium,  then  it  is  a  constant  of  the  motion  in  the  subsequent 

relaxation  motion  of  the  plasma  in  the  limit  where  the  photosphere  is 

assumed  to  be  infinitely  conducting  and  infinitely  massive.  It  is 
mathematically  allowable  to  compute  equilibria  by  specifying 
However,  if  is  increased  from  that  of  an  equilibrium,  the  coronal  plasma 
cannot  hold  it  fixed  during  relaxation  to  a  new  equilibrium.  It  is  not  a 
constant  of  the  motion  under  ideal  MHD  motion.  That  is,  a  specified  value 
of  B^  need  not  be  attainable  from  the  initial  equilibrium.  If  one  is  to 
study  quasi-static  response  to  a  slow  photospheric  motion,  slow  enough  so 
that  the  coronal  inertia  is  negligible,  then  the  footpoint  displacement, 
rather  than  B^,  must  be  specified. 

The  internal  energy  of  a  volume  of  plasma  can  be  represented  by  the 
pressure  p  or  entropy  s  (to  be  defined  below).  If  the  entropy  of  an 

initial  equilibrium  is  increased,  then  it  is  conserved  in  the  relaxation 

process  in  the  limit  where  the  plasma  is  assumed  to  be  adiabatic  on  the 

relaxation  time  scale.  That  is,  there  exists  a  limit  in  which  the  plasma 

can  be  described  by  an  equation  conserving  the  entropy  (eq.  [6dJ). 
However,  no  conservation  law  based  on  a  reasonable  energy  equation  has  been 
identified  in  which  the  pressure  is  conserved  and  which  is  suitable  for  the 

coronal  environment.  A  third  limit  which  may  be  relevant  to  the  corona  is 

the  isothermal  limit. 

In  reality,  coronal  plasmas  need  not  correspond  strictly  to  any  of  the 
limiting  cases,  in  which  case  none  of  the  physical  quantities  discussed 
would  be  specifiable  rigorously.  Specification  of  a  quantity  is 
physically  meaningful  only  in  a  suitable  limit  in  which  the  quantity  is  a 
conserved  constant  of  the  motion.  Nevertheless,  it  may  still  be  useful  to 
adopt  one  or  more  of  the  limiting  assumptions  and  explore  theoretical 
issues.  In  the  next  section,  we  discuss  the  limiting  conditions  under 
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which  the  above  quantities  are  conserved.  Additional  physical  implications 
of  specifiabili ty  are  discussed  in  Sec.  IV. 


c)  Footpoint  Displacement  and  Entropy 

By  integrating  the  z-component  of  the  field  line  equations  dx/ds  =  fi/B 
(ds  measures  arc  length  along  the  magnetic  field)  from  footpoint  1  to 
footpoint  2,  we  find 


d(.^)  =  B2(i<»)V' («<»), 


(6a) 


where 


V'(<j/) 


i: 


ds/B 


(6b) 


Here,  V(»|/)  =  Jdxdy  =  JdSpd*!// 1  =  Jdsd^/B,  where  dsp  measures  arc  length 
in  the  transverse  plane.  The  quantity  V(^)  is  proportional  to  the  area  in 
the  x-y  plane  enclosed  by  the  flux  surface  labelled  by  ^  and  V’(^)  =  dV/d^, 
so  that  V'(\|/)d»|/  is  the  volume  per  unit  length  in  the  z  direction  between 
the  flux  surfaces  labelled  and  li/  +  d\|>.  Also,  if  we  define  ♦(<!')  =  /dxdyB^ 
to  be  the  axial  flux  bounded  by  the  flux  surface  labelled  by  and  the 
photosphere,  then  we  find  from  eq.  (6a) 

=  d^  • 

During  the  process  of  increasing  d(V')  with  respect  to  time,  the  velocity  v^ 
produces  a  tangential  electric  field  =  v^By  on  S,  which  is  responsible 
for  the  change  of  ♦  in  time.  Thus,  the  footpoint  displacement  is  directly 
specifiable  through  the  boundary  conditions  and  is  a  constant  of  motion  in 
any  process  in  which  the  plasma  relaxes  to  equilibrium  on  a  time  scale 
faster  than  the  time  scale  for  the  change  in  footpoint  displacement.  We 
point  out  that  Clebsch  variables  have  also  been  used  in  specifying 
footpoint  displacement  (e.g.,  Yang,  Sturrock,  and  Antiochos  1986; 
Klimchuck,  Sturrock,  and  Yang  1988).  See  Appendix  A. 
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If  we  assume  that  the  coronal  gas  is  an  ideal  gas,  then  the  entropy 
ff'('l')d<J;  s  (da/d<J»)d<(»  of  this  volume  is  such  that 


a'(i|;)  =  Qln(pV'Y) 

where  the  contribution  YQln(d«J;)  has  been  dropped  because  only  the 
difference  in  entropy  is  significant.  Here,  y  =  5/3  is  the  adiabatic 
index.  For  convenience,  we  will  refer  to  the  following  quantity  s(<|/)  as 
"entropy"  in  this  paper  but  the  above  relationship  is  understood; 


s(4»)  =  P(»|/)(V'(<|;)1Y  . 


(6c) 


Then  s(\f>)  =  exp(  a' (i;»)/QJ . 

In  our  formulation,  the  entropy  s(\j>)  replaces  the  pressure  p(<J»)  as  the 
specified  quantity.  It  can  be  shown  that  the  entropy  is  conserved  in  MHD 
directly  without  resorting  to  thermodynamic  considerations.  Specifically, 
the  entropy  is  conserved  by  adiabatic  motions  satisfying  the  adiabatic  law 


at 


+  v7p  +  YpV'V  =  0. 


(6d) 


This  is  easily  seen  by  noting  that  p^'^Y  satisfies  the  same  continuity 
equation  as  the  density  an/at  +  7'(ny)  =  0.  Then,  the  quantity  S  a  Jp^^Y^y 
is  conserved  if  there  is  no  transport  of  plasma  across  the  flux  surfaces. 
This  condition  holds  if  the  plasma  satisfies  ideal  MHD  so  that  the  flux 
surfaces  move  the  the  fluid.  Here,  the  integral  is  over  the  volume  between 
any  two  flux  surfaces.  If  we  take  a  small  volume  6V  between  two  closely 
spaced  flux  surfaces  at  ^  and  iff  +  then  5V  is  proportional  to  V'(i^)8v». 
Therefore,  S  is  equal  to  p('|»)^'^YV'(\j()5!|;.  Since  is  conserved  for  ideal 
MHD  motions,  the  entropy  s(<J>)  a  SY/Sg/Y  is  conserved  for  an  adiabatic 
process.  Grad,  Hu,  and  Stevens  (1975)  have  developed  algorithms  to  specify 
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entropy  in  toroidal  systems.  In  the  arcade  geometry  treated  here,  an 
additional  assumption  of  negligible  heat  flux  to  the  photosphere  is 
required.  As  in  the  justification  for  specifying  footpoint  displacement, 
the  entropy  is  the  directly  specifiable  response  to  varying  the  internal 
energy  of  the  plasma  in  the  adiabatic  limit.  Note  that  pressure  is  not 
conserved  in  adiabatic  processes.  In  the  following  discussion,  we  will 
adopt  the  scenario  in  which  the  internal  energy  of  the  plasma  is  specified 
by  specifying  the  entropy  s(»(/). 

In  the  isothermal  limit,  the  temperature  remains  fixed  (perhaps  by  some 
radiative  proces)  during  the  relaxation  to  equilibrium.  An  isothermal 
plasma  obeys  a  relationship  similar  to  eq.  (6c).  If  the  mass  M(<|») 
contained  under  each  flux  surface  is  conserved  (no  photospheric  sources), 
then  M'(ij;)  =  mn( ij/)V' ( ij»)  may  be  specified,  where 

M(g»)  =  m|n(4»)^|^  . 

Mathematically,  it  is  simple  to  see  (Appendix  B)  by  comparing  eqs.  (6a) 
and  (6c)  that  p(ii/)  can  be  identified  with  and  s(i|/)  with 
(l/2)d(i</)^,  having  y  =  2  and  y  =  5/3,  respectively,  a  rather  insignificant 
difference.  For  the  isothermal  case,  s(<i/)  is  replaced  by  M'(\j>)  with  y  =  1. 

d)  Simple  Analytic  Equilibria 

We  now  illustrate  some  basic  physical  features  of  equilibrium  solutions 
of  eq.  (5)  using  force-free  examples.  A  particularly  simple  example  can  be 
given  by  the  form  B2(i|/)  =  where  Xp  is  a  constant.  Then,  the  G-S 
equation  (eq.  I5a])  reduces  to  the  Helmholtz  equation 


72^  +  X2,,,  =  0. 


(7) 


Solutions  confined  to  0  <  x  <  a  with  »|»(0,y)  =  <|»(a,y)  =  0  but  extending  to  y 
=  00  (0  <  y  <  <”)  have  the  form 
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•~k  V 

g;  =  [— Je  ,  (8) 

with  kj,  =  l(nii/a)^  -  This  equilibrium  is  the  two-dimensional, 

linear  arcade  analogue  of  the  Lundquist  (1951)  solution  (one-dimensional, 
cylindrical  geometry)  and  has  been  used  by  Heyvaerts  and  Priest  (1984)  in  a 
model  for  coronal  heating.  If  the  flux  on  the  surface  S  is  taken  to  be 
\J/(x,0)  =  sin(iix/a),  then  only  n  =  1  occurs.  In  this  case,  the  imposed  Xq 
must  satisfy  k^  >  0  or  Xq  <  X^  s  n/a.  That  is,  if  it  were  physically 
reasonable  to  specify  Xq  IB2('J')  =  »  then  one  might  expect  to  encounter 

a  violent  loss  of  equilibrium  as  Xq  is  increased  past  Xq  =  n/a. 

Using  eq.  (6a),  we  obtain 


2 

d(«(.)  =  B  («;)[ 

'’Xj  X 

(9a) 

2X 

(9b) 

where  Xq  =  a/2  is  the  point  where  ^  is  maximum  (<|/  =  ^I/q)  on  S  with  By  =  0, 
and  X  =  x(\(»)  is  the  inverse  of  =  ^(x,0).  Also,  xj^  and  X2  are  the  x 
values  where  the  flux  surface  labelled  by  ^  crosses  the  photosphere.  The 
footpoint  displacement  d  is  linear  in  x  and  the  maximum  in  magnitude  of  d 
occurs  at  X  =  0,  x  =  a  and  equals 


d 


max 


X  a 
o 


^o] 


(10a) 


Since  the  form  of  d(i^)  is  unchanged  by  Xq,  we  can  consider  d^ja^  to  be 
specified  and  Xq  to  be  computed  after  the  fact: 


(d. 


X  = 
o 


max 


(  max) 


1/2 


(10b) 
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Since  the  maximum  of  (B2)niax»  equals  consider  eq.  (9)  to 
be  a  computation  of  d(\</)  when  is  specified  and  eq.  (10b)  to  be  a 
computation  of  B^C'i/)  when  d(tj/)  is  specified.  In  the  latter  case,  no 
equilibrium  catastrophe  is  reached;  as  d,„gjj  ->  ®,  Xq  merely  increases, 
approaching  the  limiting  value  =  n/a.  This  model  is  illustrative, 
especially  since  it  also  shows  the  relationship  between  -»  ®  and 
vertical  flux  surface  expansion,  which  is  equivalent  to  -»  0  (see  eq. 
[8]).  However,  it  has  two  special  properties  that  do  not  hold  in  general; 
(i)  Xq,  or  equivalently  (B2)niax*  rnonotonicaily  increases  with  and  (ii) 
the  form  of  B2('J>)  does  not  change  as  changes  with  the  form  of  d(<|/) 
fixed. 
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III.  Linear  Arcade  Equilibria 
a)  Some  General  Properties 

In  this  section,  some  general  properties  of  force-free  configurations 
are  discussed.  First,  we  point  out  that  an  arcade  equilibrium  on  the  half 
space  -®  <  X  <  <»,  y  >  0  with  flux  and  footpoint  displacement  d(ij/) 

specified  on  the  x-axis,  must  have  a  surrounding  region  of  potential  field 
with  X((</)  =  dB^/dvJ/  =  0  or  at  least  with  and  dB^/d*//  approaching  zero 

sufficiently  fast  as  (x,  y)  -»  “.  The  footpoint  displacement  in  the  our 

model  is  such  that  d  ^  0  as  x  -»  +®,  representing  an  arcade  whose  footpoint 
shear  is  confined  to  a  finite  extent  in  x  in  the  photosphere.  For  this 
type  of  arcades,  X(ij»)  cannot  be  constant  everywhere  or  asymptotically 
approach  a  non-zero  constant.  This  can  be  seen  by  the  following 
consideration.  Here,  we  specialize  to  the  case  in  which  (J»(x,0)  is  an  even 

function  about  '  =  0,  monotonically  decreasing  for  x  >  0.  Without  loss  of 

generality,  we  may  demand  that  ^(x,0)  0  as  x  ■*  ±®.  As  a  concrete 

example,  assume  B2('J»)  =  Xq4'  +  Bq,  where  Bq  is  the  B^  field  at  infinity. 
Suppose  Xq  ^  0  so  that  dB^/dt);  ^  0  for  all  x.  Now,  suppose,  for 
contradiction,  that  Bq  it  0.  Then,  V'(ip)  becomes  infinite  for  ->  0  as  (x, 
y)  -*  ®.  This  implies  (eq.  (6aJ)  that,  if  Bq  ^  0,  then  an  infinite 
footpoint  displacement  d(<J/)  must  result  as  x  +®  and  ^  0.  Therefore, 

the  outermost  flux  surfaces  must  have  Bg  =  Bq  =  0  at  x  =  ±®  if  d(<|»)  is  to 

vanish.  Thus,  we  set  Bq  =  0.  Then,  eq.  (5a)  leads  to  +  X^tJ/  =  0.  If 
<(>(k.)  is  the  Fourier  transform  of  '(((XjO),  the  solution  takes  the  form 

00 

'l'(x»y)  =  r  dk+(k)cos(kx)e~^^^^^  (11) 

_00 

where  u(k)  =  (k^  -  X^)^^^.  Note  that 


dx4>(x,0)cos(kx)  =  2ni(>(k). 

'*—00 

Since  ij/(x,0)cos(kx)  is  even  in  x,  4'(k)  is  nonzero  in  general.  In 
particular,  ♦(k)  is  nonzero  for  k  <  Xq.  The  solution  must  then  be 
oscillatory.  This  physically  unacceptable  behavior  occurs  because  of  the 
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supposition  Xq  ^  0.  For  a  more  general  case  with  arbitrary  the 

above  argument  shows  that  X(4»)  =  must  approach  zero  as  0.  If, 

say  X(>p)  is  assumed  to  be  a  constant  X^  ^  0  in  some  region  y  >  y^*  then  the 
above  argument  applies  to  the  region  -®  <  x  <  y  >  y^,  again  with  the 

result  that  \|/  cannot  be  well  behaved  for  y  ->  ®.  This  is  a  generalization 
to  the  half-space  with  line-tying  of  the  well  known  virial  theorem  that  a 
plasma  cannot  be  confined  by  its  own  fields-  As  anotlier  example  of  virial 

theorem  application,  we  note  that  if  there  is  no  potential  field 

surrounding  the  localized  i  ^  0  region,  then  the  structure  cannot  be  force- 
free  and  the  pressure  must  be  lower  in  that  region  than  in  the  surrounding 
region.  This  has  been  shown  for  a  three-dimensional  "toroidal"  solar 

current  loop  structure  (Xue  and  Chen  1983).  Note  that  an  earlier  example 
of  a  force-free  field  without  a  surrounding  potential  field  (Low  1977)  is  a 
configuration  in  which  the  footpoint  displacement  has  an  infinite  extent  in 
the  photosphere  and  asymptotically  approaches  a  finite  maximum  at  x  =  +®. 

b)  Effects  of  Finite  Computational  Domains 

As  mentioned  earlier,  the  boundary  conditions  are  that  <|/(x,0)  is 

specified  for  0  <  x  <  a  on  S,  and  =  0  on  the  remaining  three  surfaces; 
the  two  side  walls  at  x  =  0,  x  =  a,  and  the  top  wall  at  y  =  b.  The  flux 
v(/(x,0)  is  chosen  to  be  even  about  x  =  a/2  so  that  the  maximum  value 

occurs  at  X  =  a/2  and  decreases  monotonically  to  zero  at  x  =  0  and  x  =  a. 

We  consider  a  current-carrying  plasma  confined  away  from  the  side  and  top 

walls.  Let  the  outermost  flux  surface  of  the  current-carrying  region  be  V' 
=  This  surface  intersects  the  surface  S  at  x  =  Xy  and  x  =  a-x^  where 

Xy  =  x(\^y).  The  flux  surfaces  outside  this  region  (ij/  <  'I'y)  are  then 

potential  flux  surfaces.  Based  on  the  arguments  of  the  preceding  section, 
we  specify  d(<J/)  to  be  zero  (thus  B^  =  0)  on  the  potential  flux  surfaces. 
For  =  0,  the  current  fills  the  entire  box.  Then  a  force-free 

configuration  is  supported  by  the  walls  and  the  form  of  B2('l')  for  specified 
d(<f/)  will  be  seriously  modified.  For  larger  a  force-free  structure  is 
supported  by  the  surrounding  potential  flux,  which  in  turn  is  supported  by 
the  walls.  If  is  sufficiently  large,  then  the  interface  surface  g;  =  (1;^ 
is  sufficiently  far  from  the  walls  and  its  shape  is  determined  by  the 
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actual  force  balance  requirement.  In  this  section,  we  will  discuss  the 
generic  effects  of  the  finiteness  of  the  computational  domain  on  the 
solution.  Hereafter,  we  use  dimensionless  variables  with  x,y,  d(<|/)  scaled 
to  a  (0  <  X  <  1).  The  flux  function  is  scaled  to  its  maximum  is 
scaled  to  'J/g/a,  and  p  is  scaled  to 

We  expect  the  general  form  of  d(<j»)  to  be  linear  in  x  near  the  neutral 
line  at  X  =  a/2  [as  it  is  in  eq.  (9b))  and  vanish  as  the  potential  region 
is  approached.  Clearly,  d(\|/)  must  vanish  at  x  =  a/2  because  the  field  line 
has  zero  length  there.  Furthermore,  if  the  magnetic  field  is  regular  in 
that  vicinity,  then  d  is  linear  in  x  unless  is  finite  and  both 
components  and  By  vanish  at  x  =  a/2,  y  =  0.  Otherwise,  an  0-point  rises 
above  the  boundary  surface  S.  This  is  inconsistent  with  specifying  d(\J/); 
if  an  0-point  forms,  then  d  must  have  a  step  function  discontinuity. 

As  a  first  example,  we  have  computed  force-free  equilibria  with  d(x)  of 
the  form 


d(x)  =  d^ 


tanh 


r  -  \>  1 

l»(a'3  -  x^)J 


X  <  X  <  a/2 
V 

('I'  >  ’I'y). 


(12a) 


=  0 


X  <  X 

V 

('<'  <  'I'y) 


(12b) 


and  anti-symmetric  about  x  =  a/2.  The  G-S  equation  is  solved  iteratively 
for  specified  d(i/').  The  computational  method  is  described  in  Appendix  B. 
In  Figure  2a,  we  show  the  flux  surfaces  for  a  representative  equilibrium. 
We  have  used  \|;(x,0)  =  \|/Qsin( iix/a)  with  =  1,  the  height  of  the 

integration  box  b/a  =  2,  and  =  0.195,  for  which  x^  =  0.0625.  The 

maximum  footpoint  displacement  is  specified  to  be  =  1.85  and  we  have 

chosen  W  =  0.1.  The  coefficient  dg  is  adjusted  so  that  the  specified  value 
of  d,„gx  is  obtained.  The  short-dashed  contour  corresponds  to  \j;  =  where 

=  0  and  d  =  0.  All  the  current  is  contained  inside  (\(/  >  <J/y)  the  g;  = 
surface.  Thus,  the  outer  boundary  surface  of  the  i  ^  0  region  is  removed 
from  the  computational  box  by  a  moderate  amount.  The  long-dashed  contour 
corresponds  to  the  flux  surface  defined  to  be  the  flux  surface  such 

that  d(<)«y)  =  dnigjj  /2.  (For  this  example,  we  have  d(\J;y)  =  0.92,  Xy  =  0.08, 
and  ij>y  =  0.25.)  The  bulk  of  the  current  is  actually  enclosed  by  this 
surface,  farther  away  from  the  walls  than  the  surface.  Figure  2(b) 
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shows  the  footpoint  displacement  profile.  (This  figure  shows  d  only  for  0 
<  X  <  0.5.  The  overall  sign  is  reversed  from  that  of  eq.  {12al  for 
convenience.)  The  above  form  has  the  general  feature  that  the  anti¬ 
symmetric  ^'ootpoint  displacement  is  zero  for  0  <  x  <  x^  with  most  of  the 
footpoint  displacement  shear  concentrated  in  x^  <  x  <  x,^  where  =  x^  + 
U(a/2  -  x^,)  is  where  the  maximum  footpoint  displacement  d^ax  =  ^^^m^ 

occurs.  In  Xp,  i  X  ^  (a  -  x^),  d(x)  is  essentially  linear  so  that  field 

lines  on  different  flux  surfaces  are  parallel  and  unsheared.  For  larger 
values  of  W,  the  footpoint  displacement  shear  is  less  concentrated. 

In  Fig.  2(c),  we  have  plotted  B2('J')  versus  ij;  for  0  <  ij/  <  1 .  (Recall 
that  on  the  surface  S,  4/  =  0  for  x  =  0  and  x  =  1,  and  4;(x=0.5,  y=0)  =  4/^  = 

1.)  Ue  see  that  82(4^)  is  linear  in  4;  near  4/q  and  deviates  from  linearity 
in  the  region  of  maximum  shear  in  the  footpoint  displacement.  In  Fig. 
2(d),  we  show  X(4;)  =  as  a  function  of  4/.  As  expected,  X(4;)  has  its 

maximum  variation  in  the  region  of  maximum  shear  in  d(4/). 

In  Fig.  3,  the  quantity  (B2)max  “  B2(x=0. 5,y=0)  is  plotted  as  a 

function  of  for  a  number  of  model  equilibria.  Curve  1  is  a  plot  of 

eq.  (10b),  a  linear  force-free  configuration  confined  horizontally  to  0  i  x 
i  1  but  extending  vertically  to  y  =  ®.  For  this  system,  (B2)niax  Xq) 

is  linear  in  for  small  d^^j^  and  asymptotically  approaches  n/a  as  d^jgj^ 

®.  For  small  d^gj^,  V'(4/)  is,  to  lowest  order,  the  value  for  the 
potential  field  (d^jgj^  =  0)  and  (B2)niax  proportional  to  d^gj^  (see  eq. 

[6a]).  For  d^jax  ”»  other  hand,  the  increase  in  d^jgj.  causes  an 

expansion  of  flux  surfaces,  i.e.  an  increase  in  V'(4;),  with  (B2)niax 
approaching  a  constant.  Curve  2  of  Fig.  3  is  for  the  solution  of  eq.  (7) 
for  a  finite  height  b,  where  the  boundary  condition  is  4/(x,  y=b)  =  0.  The 
equilibrium  can  be  described  analytically  by 

4'(x,y)  =  sin^^jsinh  k.^(b-y)/sinh  k^b,  (13) 

where  kj  =  [(rt/a)^  -  Xq]^^^  for  n/a  >  Xq.  For  X^  >  n/a,  eq.  (13)  is 
replaced  by 

4/(x,y)  =  sin|^^jsin  kj(b  -  y)/sin  kj^b,  (14) 

where  kj  =  [Xq  -  (n/a)^!^'^^.  For  this  equilibrium  of  finite  height,  the 
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functional  form  of  d(i</)  (not  merely  the  amplitude  d^gjj)  changes  with  so 
that,  unlike  the  infinite  height  equilibrium  given  by  eq.  (8)  or  eq.  (10b) 
(curve  1),  specifying  Xq  or  (B2)[nax  essentially  equivalent  to 
specifying  For  the  infinite  height  case,  eq.  (9b),  the  profile  of  d 
remains  linear  in  x  and  that  of  linear  in  g/  for  all  values  of  Xq,  or 
equivalently,  for  all  values  of  d^gj^.  Nevertheless,  the  curves  1  and  2  are 
similar.  Curve  2,  however,  shows  that  (Bz)niax  increase  beyond 
the  limit  Xq  -*  it/a,  in  contrast  with  curve  1.  For  curve  2,  the  finite 
height  (b)  of  the  box  prevents  the  flux  surfaces  from  expanding  freely  as 
dmax  ”•  case,  V'(\i/)  cannot  increase  without  bounds  and  dj^gj^  can 
increase  only  through  an  increase  in  (^z^max  (see  eq.  [6aJ).  However,  for 
this  equilibrium,  the  flux  function  ip  takes  the  form  of  eq.  (14)  for  X  > 
n/4  and  an  0-point  appears  at  x  =  a/2  above  the  surface  S  (y  >  0)  if  kjb  > 
n/2,  or  Xq  >  (n/a)(l  +  a^/4b^)^^^.  For  b/a  -  2,  the  parameter  used  in 
curves  2  and  3,  this  gives  X  =  3.24.  We  have  not  attempted  to  extend  d^g^ 
past  this  point.  Curve  3  is  the  equilibrium  for  a  footpoint  displacement 
profile  given  by  eq,  (12a)  with  =  0.195  and  W  =  0.1,  the  equilibrium 
shown  in  Fig.  2a.  This  equilibrium  is  numerically  computed  in  the  same 
finite  box  as  for  curve  2,  As  in  curve  2,  increases  past  n/a  as 
dmax  increases,  albeit  slowly,  because  of  the  finite  size  of  the 
computational  domain.  Presumably,  (^z^max  continues  to  increase,  at  a 
faster  rate  with  respect  to  d^^gjj,  as  d^g^  increases  further.  However,  the 
non-zero  current  region  (inside  the  short-dashed  flux  surface  in  Fig.  2a) 
is  detached  from  the  side  walls  (and  the  region  <  'ky  has  less 
current  than  in  Fig.  2a).  Thus  there  is  relatively  more  room  to  expand 
than  in  the  configuration  described  by  curve  2.  This  flexibility  for  the 
flux  surfaces  to  expand  against  the  potential  flux  surfaces  is  responsible 
for  the  smaller  values  of  (B^lmax  large  d^g^.  The  difference  in  the 
slopes  of  curve  3  and  curve  1  (or  curve  2)  is  not  physically  significant, 
but  is  related  to  the  fact  that  d^gj^  itself  has  a  different  meaning  for  the 
two  very  different  profiles  of  d('k),  given  by  eq.  (9b)  and  eq.  (12).  Curve 
3  shows  that  the  numerically  computed  equilibrium  solutions  whose  current- 
carrying  regions  are  only  slightly  removed  from  the  walls  are  similar  to 
the  analytic  solutions  in  which  the  current-carrying  plasmas  fills  the 
finite  box.  The  main  effect  of  the  walls  is  to  cause  (B2)n,ajj  to  increase 
past  the  asymptotic  limit  Xq  =  n/a  of  eq.  (10b)  as  d^g^  ->  <”. 
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c)  Equilibrium  Solutions  with  Horizontal  Expansion 


The  flux  surfaces  in  Fig.  2a  have  expanded  upward  to  a  large  degree 
with  fairly  large  d^gj^  =  1.85.  In  order  to  investigate  equilibria  with 
the  walls  farther  away,  allowing  more  room  to  expand  horizontally  as  well 
as  vertically  in  response  to  increasing  we  consider  a  more  general 
class  of  equilibria  with 


<(/(x,  0)  =  sin''(iix/a)  (15) 

with  V  >  1.  For  purposes  of  illustration  we  again  consider  solution  of  eq. 
(7),  satisfying  the  boundary  condition  (19).  Then, 
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(16) 


where  k^  =  ((nn/a)^  -  If  k^  <  0,  hyperbolic  sine  is  replaced  by 
sine.  Here,  are  the  Fourier  sine  coefficients  of  eq.  (15).  In  Fig.  4, 
(Bz)max’  which  is  equal  to  Xq  (from  =  Xq\J/  and  =  1),  is  shown  as  a 
function  of  for  \)  »  7.  Ue  see  that  (B2)p,3j^  increases  past  n/a  for 
dmax  ^  2  with  a  much  smaller  slope  than  for  dp,gjj  ^  2.  Also  plotted  in 
Fig.  4  is  the  total  current  in  the  z-direction,  I^  =  Jdxdyj^  = 
JdxdyB2( +)dB2/dv|/.  This  component  of  current  is  seen  to  increase 
quadratically  for  small  dpjg^  (I2  ®  X^)  and  to  continue  to  increase  for 
larger  dp,gj^.  This  behavior  appears  to  be  general  for  straight  two- 
dimensional  arcade  equilibria,  as  will  be  shown  with  more  general  examples. 
For  the  equilibrium  given  by  eq.  (8)  with  n  =  1,  I2  is  proportional  to 
x2/ki  =  x2/(ii2/a2  -  x2)l/2. 

In  Fig  5(a),  we  show  the  flux  surfaces  for  an  equilibrium  with  v  =  7, 


2  r  ^v^  1 

d(x)  =  d^(x  -  a/2)  tanh  [v(a72  -  x  ^  ^  ^  (17a) 

=0  X  <  x^  (17b) 

(and  again  antisymmetric  about  x=a/2)  having  x^  =  0.31  =  0.25),  dp,gjj  = 


20 


0.36,  and  W  =  0.40.  For  this  example,  =  0.40  which  intersects  the 

surface  S  at  x  ^  0.34.  The  change  in  form  between  eqs.  (12)  and  (17) 

(providing  continuous  derivative  d'(x)  at  x^],  and  the  larger  value  of  W, 

serve  to  make  the  solutions  less  sensitive  to  behavior  near  =  i(/y.  The 

corresponding  profiles  of  d((J/),  B2('|/),  and  X(<(/)  are  shown  in  Figs.  5(b)- 
5(d).  The  equilibria  shown  in  Figs.  6  and  7  have  larger  footpoint 
displacements  d^^j^  =  0.61  and  1.35,  respectively.  It  is  easily  seen  that 
the  flux  surfaces  expand  upward  and  outward  as  d^gj^  increases. 

In  the  example  of  eqs.  (8),  (13)  and  (16),  and  in  the  numerical  results 
shown  in  Figs.  3  and  4,  (^z^max  continues  to  increase  with  dn,^^ 
approaching  a  limiting  value  or  increasing  depending  upon  whether  the 
upward  expansion  of  flux  surfaces,  and  therefore  further  increase  of  V'(iJ/), 
is  limited  by  a  wall  at  y  =  b.  For  the  examples  shown  in  Figs.  5-7,  the 
flux  surfaces  are  allowed  to  expand  horizontally  as  well  as  vertically. 
The  effect  of  allowing  horizontal  spreading  of  the  flux  surfaces  manifests 
itself  in  the  graph  on  (B2)n,ax  versus  d,^g^  in  Fig.  8.  For  equilibria  of 
the  series  (B)  shown  in  Figs.  5-7,  for  which  x^  =  0.31  and  W  =  0.40,  there 
is  a  peak  in  (B2)max  ''max  =  0.35.  For  another  related  series  (A)  of 

equilibria,  having  x^  =  0.17  (\J;^  =  9  x  10“^)  and  W  =  0.45,  the  peak  is 

broader  and  is  located  near  d^g^  =  0.5.  The  observed  decrease  in  (B2)niax 

is  related  to  horizontal  spreading  for  the  following  reason:  if  no 
horizontal  spreading  occurs,  as  in  eqs.  (8),  (13)  and  (16),  B2  increases 
very  slowly  for  large  d^jg^?  and  in  the  case  of  eq.  (8),  has  an  asymptote 
(Bz)n,ax  have  discussed,  this  behavior  is  due  to  the  fact  that 

vertical  spreading  in  flux  surfaces,  i.e.  the  increase  in  the  V'(iJ>)  term  in 
eq.  (6a),  can  provide  the  increase  in  d^g^  with  little  change  in  (B2)niax‘ 
When  horizontal  spreading  is  allowed  by  having  the  photospheric  flux 
concentrated  as  in  eq.  (15),  V'(ij;)  can  increase  at  a  rate  faster  than  d,„gx» 
giving  a  decrease  in  (B2)niax‘  case  of  Fig.  4  described  by  eqs. 

(15)  and  (16),  this  effect  cannot  occur.  Mathematically,  this  is  simply 
traced  to  the  fact  that  if  Xq  =  (B2)tnax  prescribed  in  eq.  (7),  then  the 
solutions  are  unique,  and  therefore  only  one  value  of  dp,gx  is  possible  for 
each  Xq.  Physically,  this  appears  to  arise  from  the  fact  that  for 
the  outer  surfaces  spread  out  horizontally  even  for  very  small  dj^g^.  when 
(Bz)max  increasing  rapidly.  Thus,  it  appears  that,  in  order  for  (B2)n)ax 
to  attain  a  peak  and  decrease  thereafter,  a  certain  amount  of  potential 
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flux,  or  at  least  a  profile  having  X(4<)  0  as  the  wall  is  approached,  is 

required  to  separate  the  current-carrying  region  and  the  fixed  walls  so 

that  the  flux  surfaces  can  spread  out  against  a  "soft"  wall  when  the 
displacement  is  large.  Nevertheless,  for  even  larger  (B2)n)ax 

again  begin  to  increase,  when  the  expansion  of  flux  surfaces  is  arrested  by 

the  walls.  This  effect  may  be  beginning  to  occur  for  d^jgj^  ^  1  in  Fig.  8. 

However,  this  effect  is  due  to  the  unphysical  boundaries.  We  expect  that 
(^z^max  continue  to  decrease  in  the  absence  of  such  boundaries. 

In  our  iteration  scheme  to  solve  the  G-S  equation  for  specified  d^g^ 
(Appendix  B),  the  outer  iteration  loop  described  by  eq.  (12)  converges  very 
rapidly  near  the  peak  in  (Bj,)^gx  beyond,  because  the  solutions  have  a 

relatively  weak  dependence  upon  d^gj^.  However,  the  inner  loop  given  by  eq. 
(Bl)  does  not  noverge  for  solution  near  and  past  this  in  (B2)|nax'  This  is 
because  B^i'k)  is  specified  in  the  inner  loop  as  an  intermediate  step  and 
the  peak  corresponds  to  a  bifurcation  in  this  case.  Although  such 
bifurcations  are  not  physical  as  discussed  in  Secs.  I  and  Ilb,  they  have 
mathematical  and  numerical  significance  for  our  iteration  scheme.  In 
general,  an  iteration  as  in  eq.  (Bl)  is  expected  to  converge  for  one  class 
of  solutions  (with  smaller  d^gj^)  and  diverge  for  the  other  class  (with 

larger  d^gj^).  We  overcome  the  divergence  problem  by  renormalizing  B^C'j')  in 
the  inner  loop  iterations  (that  is,  changing  (B2)niax  no  change  in  the 

profile  B2('J')/(B2)n,axl  order  to  conserve  the  total  z-current  = 

JdxdyB2dB2/d»J/.  This  is  a  simple  operation  because  I^  scales  as  (B|)|ngx' 
From  Fig.  4,  we  see  that  I^  increases  much  more  rapidly  with  d^gj^  than 
(B^lmax  solution  given  by  eq.  (16).  In  fact,  it  is  found  that  I2 

increases  for  all  solutions  even  when  (Bz)niax  decreases.  With  this 
prescription  in  which  I^  rather  than  (Bz)max  held  fixed  during  the 

iteration  process,  there  is  no  bifurcation.  When  this  method  is  employed, 
the  inner  loop  iterations  are  always  found  to  converge. 

For  d^gj^  =  0.36  (Figs.  5a  -  5d),  the  region  of  strong  shear  in 
footpoint  displacement  is  0.31  <  x  <  0.34,  or  equivalently  0.26  <  <  0.40. 
In  this  region,  Bzit)  [Fig  5(c)j  deviates  significantly  from  linearity. 
Note  that  X(<j/)  =  dB^/dij/  is  not  constant  (Fig.  5d),  decreasing  to  zero  as  <1^ 
decreases.  This  behavior  is  consistent  with  the  general  properties 
discussed  in  Sec.  Ilia.  Closer  to  the  center,  x  >  0.34  and  >  0.40,  Bz  is 
fairly  linear  in  ><;,  although  X  shows  a  slight  peaking  at  the  center.  Since 
i  =  X(ij;)B  for  these  equilibria,  this  indicates  a  peaking  of  the  current 
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density  relative  to  the  solutions  of  the  form  eqs.  (8),  (13)  and  (16).  In 
Fig.  6,  for  a  larger  footpoint  displacement  with  d^^jj  =0.61,  this  tendency 
appears  to  be  somewhat  more  pronounced.  Although  not  evident  in  Fig.  6(b), 
Figure  6c  shows  a  sharp  drop-off  in  X  in  the  shear  region  0.26  <  ^  <  0.40, 
a  further  peaking  near  the  center  (t|/  >  0.8)  and  a  plateau  in  between.  For 

a  still  larger  footpoint  displacement  dj^gj^  =  1.35  (Fig.  7c),  this  tendency 
is  even  more  pronounced,  with  X(»|/)  having  a  rather  narrow  peak  for  \)/  >  0.9, 
sharp  drop-off  in  the  shear  region  0.26  <  \|/  <  0.40,  and  a  valley  around  = 
0.7.  The  current  density  ~  quite  peaked;  its  value  at 

the  center  x  =  a/2,  y  =  0  (»|/  =  1)  is  65  (in  normalized  units),  compared  to 

a  considerably  lower  value  of  10  at  <//  =  0.7.  That  is,  in  a  sense, 

solutions  with  increasing  d^jgjj  exhibit  a  considerable  amount  of  peaking  in 
the  current  density.  Nevertheless,  we  will  show  below  that  the  peaking  in 
at  the  center  with  increasing  d„,gjj  does  not  lead  to  formation  of  current 
sheets.  Moreover,  the  current  density  perpendicular  to  the  photosphere  jy 
=  XBy  =  -X3\;//3x  is  proportional  to  sin^( nx/a)cos( iix/a)X(\^) .  Because  By  is 
small  near  x  =  a/2  (<^  =  1),  the  peak  in  X  also  does  not  correspond  to  a 

large  peaking  of  jy. 

In  Fig.  9,  X^ax>  maximum  value  of  X(\|/),  is  shown  as  a  function  of 
dmax*  This  maximum  occurs  at  i)/  =  for  all  cases.  Recall  that  X^ax  =  \) 
is  identical  to  (B^lmax  solutions  eqs.  (8),  (13),  and  (16).  For  the 

equilibria  of  Fig.  9,  however,  there  are  qualitative  differences.  Note 
that  X^ax  begins  to  flatten  for  increasing  d^jgj^,  at  d„,gjj  =  0.2  for  series  B 
and  at  d^jg^  =  0.4  for  series  A  while  (Bz)max  peaks  at  x  =  0.35  and  x  = 
0.5,  respectively.  However,  in  both  cases  X^ax  continues  to  increase, 
unlike  (B^lmax'  seems  likely  that  X^g^  may  have  a  peak  with  respect  to 
dmax  larger  values  of  x^  and/or  W.  However,  for  the  present 

parameters,  the  proximity  of  the  wall  causes  the  slow  increase  in  X^ax 
dmax  ^  O.A.  This  behavior  shows  further  the  lack  of  any  physical  meaning 
to  the  "bifurcations"  in  Fig.  8. 

We  now  consider  the  peaking  in  X((|>)  and  ^^('l')  in  more  detail.  In  Fig. 
10,  we  show  the  maximum  current  density  as  a  function  of  d,„gj^.  For 

solutions  of  eq.  (7)  (for  example,  eqs.  [8],  [13),  and  [16)),  j^  is 
proportional  to  Uz  =  see  eq.(7)).  We  see  that  for  series  A  and  B, 

jj,  and  X  are  seen  to  have  similar  behavior;  increases  rapidly  for  d^jg^  i 
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0.35  (series  B)  and  for  d^jgjj  ^  0.50  (series  A),  the  same  regions  for  rapid 
increase  in  Figs.  8  and  9.  Also,  continues  to  increase,  but  much  more 
slowly,  for  larger  values  of  The  slow  increase  of  jj,  for  d,„gjj  i  0.5 
appears  to  be  due  to  the  proximity  of  the  wall.  Ue  also  expect  that 
^jz^max  '"ight  decrease  with  d^j^^  ^^r  larger  x^  and/or  W.  At  the  same  time, 
the  profile  of  =  XB2  continues  to  narrow  so  that  the  total  current  in 
the  peak,  region  decreases.  In  Fig.  11,  we  show  the  total  z-current  is 
shown  for  series  A  and  B.  It  is  observed  to  increase  quadratically  for 
*^max  ^  O'l  snd  monotonically  for  all  (as  noted  above.)  However,  I^ 
does  not  increase  as  rapidly  for  large  d^jgj.  as  in  Fig.  4.  The  important 
points  to  note  are  (i)  that  the  peak  current  density  is  always  found  to 
increase  at  a  very  slow  rate  with  respect  to  giving  a  decrease  in  the 
total  current  in  the  peak  region.  Concurrently,  (ii)  the  total  current 
I2  increases  rapidly  with  so  that  the  fraction  of  the  total  current 
I2  carried  by  the  peak  decreases  very  rapidly  with 

We  point  out  that  the  single-valuedness  of  the  graphs  in  Figs.  8-11  is 
due  to  the  fact  that  a  single  solution  is  found  when  d(ij;)  is  specified.  We 
take  this  as  empirical  evidence  that  these  solutions  are  unique. 
Therefore,  any  other  physical  quantity  such  as  the  total  magnetic  energy 
will  also  be  single-valued.  It  has  also  been  found  that  the  total  magnetic 
energy  of  a  linear  dipole  field  can  be  well  fitted  with  a  logarithmic 
dependence  on  footpoint  shear  (Klimchuk,  Sturrock,  and  Yang  1988). 

As  we  have  discussed  in  Sec.  lie  and  Appendix  B,  finite  3  equilibria 
with  entropy  prescribed  can  be  computed  in  essentially  the  same  manner.  In 
Fig.  12  we  show  results  with  dp^g^  =  0  (i.e.  =  0)  as  a  function  of  the 
maximum  entropy  As  suggested  by  the  analogy  between  eqs.  (B4)  and 
(B5),  the  form  of  s(\J/)  is  chosen  to  be  ld(\J/)J^^^,  where  d(ij;)  is  the 
footpoint  displacement  form  factor  of  eq.  (17).  Fig.  12a,  we  show  the 
dependence  of  p^ax  versus  for  Xy  =  0.31  and  W  =  0.40,  the  same  choices 
as  in  curve  B  in  Fig.  8.  We  have  used  zero  footpoint  displacement  since 
the  results  of  Zwingmann  (1987)  show  that  this  is  the  case  with  the 
strongest  degree  of  bifurcation  when  pressure  is  specified.  We  see  that 
Pmax  3  well-behaved,  single-valued  function  of  For  very  small 
s^gj^,  increasing  s^gj^,  i.e.,  the  internal  energy  content,  predominantly 
leads  to  a  linear  increase  in  p^ax*  This  curve  has  a  broad  maximum  around 
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Smax  ~  increasing  in  this  region  leads  to  expansion  of  the 
structure  in  such  a  way  that  P^ax  remains  nearly  unchanged.  For  i 
0.4,  increasing  is  predominantly  accomplished  by  expansion  of  the 
structure  and  the  pressure  in  fact  decreases.  This  behavior  is  similar  to 
that  of  (Bjrlnjgjj  versus  d^gj^  (see  Fig.  8).  The  flux  surfaces  are  similar 
to  those  shown  in  Figs.  5a,  6a,  7a.  The  overall  behavior  of  equilibria 
with  respect  to  variation  of  s^g^  is  generally  similar  to  that  with  respect 
to  dj^g^  with  no  remarkable  differences.  We  note  that  Fig.  7  of  Zwingmann 
(1987),  which  corresponds  to  a  solution  on  the  upper  branch  of  the  W 
(energy)  versus  Xp  (a  measure  of  pressure)  curve,  shows  apparent  formation 
of  an  0-point.  Such  an  equilibrium  cannot  evolve  from  configurations 
without  0-points  via  ideal  MHD  motions.  Mathematically,  such  0-points  will 
not  occur  when  d(v(;)  or  s(vj/)  is  prescribed  because  these  quantities  are 
specified  only  on  the  range  0  <  tp  <  and  the  occurence  of  an  0-point 
indicates  the  existence  of  surfaces  with  >  \)/q,  the  maximum  flux  value 
specified  on  S.  The  occurence  of  0-points  when  B2(t)  is  specified  has  also 
been  found  by  Birn,  Goldstein,  and  Schindler  (1978). 

In  Fig.  12b,  we  show  the  total  axial  current  and  the  peak  current 
density  j2(x=a/2,  y=0).  As  the  entropy  increases,  the  current  density  jg  = 
dp/di|/  becomes  peaked  at  the  center,  as  in  the  case  of  increasing  the 
footpoint  displacement.  However,  the  total  current  again  increases  at  a 
faster  rate  with  respect  to  s^gj^  than  the  current  in  the  peak  region  (x  = 
a/2),  so  that  the  fraction  of  the  current  associated  with  the  peak 
decreases  relative  to  as  s^^g^  increases.  No  remarkable  differences 
appear  to  exist  between  specifying  d^g^  and  specifying  s^gj^,  consistent 
with  the  mathematical  similarity  apparent  in  the  formulations  in  eqs.  (6a) 
and  (6c).  Note  that  p^ax  well-behaved,  single-valued  functions 
of  s^gjj.  Thus,  the  total  energy,  which  is  the  sum  of  the  magnetic  energy 
and  internal  energy  (“:p)  integrated  over  the  volume,  is  single-valued, 
exhibiting  no  bifurcations. 
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IV.  Discussion 


The  basic  physical  mechanism  behind  the  expansion  of  the  flux  surfaces 
'^max  ^max  increased  is  relatively  simple.  The  force-free 

equilibrium  condition  i  x  S  =  0  can  be  written  as  (jpBg.  -  jz®p)  =  where 
jp  is  the  transverse  (x,y)  component.  For  a  force  free  equilibrium,  jpB^  ~ 
dB|/d\j;  is  an  upward  force  if  B^  monotonically  increases  in  <</.  In  an 
imaginary  two-step  process  (Sec.  Ilb),  if  the  footpoint  displacement  is 
increased,  at  first  with  no  change  in  the  flux  surfaces  [no  change  in 
Bp  or  jzl>  B^  must  increase.  However,  the  plasma  will  be  out  of 
equilibrium  and  will  expand  upward  under  jpB^  until  force  balance  is 
achieved  again.  This  expansion  will  occur  with  dp,^^.  conserved  in  the  limit 
of  infinitely  conducting  photosphere,  which  means  that  the  flux  between 
any  two  flux  surfaces  will  be  conserved.  This  expansion  process  will  then 
necessarily  decrease  In  the  new  equilibrium,  B^  may  be  greater  than  or 

less  than  that  of  the  initial  equilibrium,  depending  on  whether  the  initial 
increase  of  B^  dominates  the  subsequent  decrease  of  B^.  There  is  a  simple 
circuit  analogy  to  this  process.  As  we  have  seen  in  the  discussion 
following  eq.  (6),  d('J»)  is  proportional  to  the  axial  flux  between  two 
nearby  flux  surfaces.  Also,  B2('l')  is  essentially  the  total  transverse 
current  between  the  flux  surface  labelled  by  and  that  labelled  by  (see 
the  discussion  after  eq.  (3J).  By  eq.  (6a),  V'(i<^)  (which  is  a  property 
only  of  the  geometry  of  the  flux  surfaces)  can  therefore  be  identified  as 
the  inductance  relating  transverse  current  to  axial  flux.  In  the  first 
step  of  the  imaginary  two-step  process  above,  when  the  photospheric  motion 
takes  place  (d^gjj  is  increased)  with  the  geometry  held  fixed,  we  have 
Ldl/dt  =  Vq.  The  voltage  Vg  is  proportional  to  the  tangential  electric 
field  on  the  photospheric  surface  Ej^  =  v^By.  In  the  second  step,  the 
geometry  and  current  change,  conserving  flux  d(LI)/dt  =  0.  For  small  dp^gj^, 
the  process  is  completed  primarily  by  increasing  the  current  with  little 
change  in  inductance;  for  large  clp,gj^,  the  inductance  increases  by  an  amount 
sufficiently  large  that  the  current  must  decrease. 

A  similar  argument  applies  to  the  pressure.  In  thermal  equilibrium, 
the  thermal  input  and  radiative  loss  between  any  two  flux  surfaces,  say, 
and  ij/2’  balanced.  If,  for  example,  the  heating  rate  increases,  the 

pressure  "initially"  (i.e.,  before  relaxation  to  equilibrium)  increases. 
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Then,  in  approaching  equilibrium,  the  flux  surfaces  and  ^2  will  both 
expand  under  7p  in  such  a  way  as  to  decrease  the  pressure,  by  adiabatic 
decompression  (assuming  that  MHD  time  scales  are  shorter  than  radiative  or 
thermal  time  scales).  During  this  process,  if  expansion  dominates,  then 
the  new  equilibrium  pressure  will  be  lower  than  that  of  the  initial 
equilibrium.  The  imaginary  two-step  process  described  in  Sec.  Ilb  and 
above  is  helpful  not  only  for  conceptualization  but  for  predicting  the 
response  of  the  coronal  plasma  to  a  source  of  heat  that  occurs  on  a  faster 
time  scale  than  the  nominal  heating  and  radiative  loss  mechanisms. 

The  scenario  for  which  entropy  rather  than  pressure  should  be  specified 
assumes  that  the  photosphere  does  not  act  as  a  source  of  material  or 
energy.  This  can  achieved  if  thermal  conduction  along  the  field  lines 
between  the  corona  and  the  photospheric  or  subphotospheric  regions  is 
negligible  on  the  MHD  time  scale  for  relaxation  to  equilibrium,  or 
equivalently,  if  the  magnetic  field  in  the  arcade  is  large  enough  and  the 
coronal  density  small  enough  that  the  Alfven  time  scale  is  much  shorter 
than  the  time  scale  for  parallel  thermal  conduction  in  the  photosphere. 
This  argument  also  assumes  that  radiation  is  negligible  on  the  Alfven  time 
scale.  On  the  other  hand,  if  the  radiation  is  dominant,  and  the  plasma  is 
isothermal,  the  mathematical  formulation  is  essentially  identical  to  the 
adiabatic  case  (but  with  adiabatic  index  y  =  !•)  In  our  computation,  we 
have  also  assumed  that  the  arcade  length  scales  are  short  compared  with  the 
coronal  gravitational  scale  height  2kT^,/mg,  so  that  p  =  p(\J/)  in  the  corona. 
This  is  only  for  convenience,  and  is  not  critical,  especially  if  p  =  2p/B^ 
is  small.  (The  photospheric  scale  height  2kTp/mg,  on  the  other  hand,  is 
assumed  to  be  small.  This  is  important  in  order  to  be  able  to  specify 
footpoint  displacement.  This  is  consistent  with  Tp  being  much  less  than 
Tj,.)  On  the  other  hand,  suppose  one  attempts  a  justification  for 
specifying  coronal  pressure  based  on  the  dominance  of  photospheric 
processes  by  specifying  the  photospheric  pressure  p(\|/,0).  Accordingly, 
assume  that  an  increased  photospheric  pressure  is  imposed  with  no  changes 
in  the  flux  surfaces,  as  in  the  imagined  two-step  process  discussed  above. 
At  this  step  in  the  process,  the  increased  coronal  pressure  is  given  by  eq. 
(4b).  Then,  the  pressure  p('l',y)  will  subsequently  decrease  as  the  system 
relaxes  to  equilibrium.  There  does  not  seem  to  be  a  reasonable 
conservation  law,  one  based  upon  an  energy  equation  for  the  plasma,  in 


which  the  pressure  in  the  corona  p(<(<,y)  given  by  eq.  (4b)  can  be  specified. 
Using  this  interpretation,  we  conclude  that  bifurcations  with  respect  to 
prescribing  pressure  do  not  correspond  to  a  physical  loss  of  equilibrium. 

The  effect  of  expansion  of  the  flux  surfaces  and  the  effect  upon  the 
curves  in  Figs.  8  and  12a  can  be  understood  better  by  the  following  scaling 
argument.  First,  consider  the  force  free  case  s(\f/)  =  p(<p)  =  0.  For 
arcades  that  are  free  to  expand  upward  but  are  constrained  in  the 
horizontal  direction  (e.g.,  eqs.  (81  and  [13];  also  Fig.  2a),  the  Grad- 
Shafranov  equation  (5a)  implies  in  the  large  d^^j^  limit  that  ~  '('/a, 
independent  of  the  height  h  of  the  flux  surfaces.  From  eq.  (6a),  we  have  d 
=  B^V'  so  that,  for  large  h,  V'  ~  ah/\l/,  d  ~  (\(//a)(ah/\p)  and  h  ~  d.  On  the 
other  hand,  if  the  plasma  is  free  to  expand  in  both  directions,  its  width 
will  be  of  the  same  order  as  its  height  h  for  large  h  and  the  G-S  equation 
gives  B^  ~  'J/Zh.  Using  V'  ~  we  again  obtain  d  ~  (t;//h)(h^/ii»)  ~  h. 

This  in  turn  implies  B^  ~  ^/d  for  large  d.  This  exact  scaling  is  not 

observed  in  the  results  of  Fig.  8  for  two  reasons.  First,  the  effect  of 
the  proximity  of  the  wall  is  probably  non-negligible  for  dp,gjj  ^  1.0  in  Fig. 
8.  Second,  such  a  scaling  argument  cannot  include  the  effects  of  profile 
changes,  e.g.,  the  peaking  in  the  X  profile  seen  in  Fig.  7c.  Nevertheless, 
this  scaling  does  illustrate  the  general  features  of  the  results  observed. 

A  similar  scaling  law  can  be  obtained  if  entropy  is  specified  with  d(»(') 
=  B2('(')  =  0.  If  only  vertical  expansion  is  allowed,  then  p  ~  (from 

the  G-S  equation)  and,  using  s  =  pV'T,  we  obtain  s  ~  (\(^/a^)(ah/\j/)'Y  for 
large  h,  or  s  ~  hT.  For  two-dimensional  expansion,  we  obtain  p  ~  and 

s  ~  (<p^/h^)(h^/\;/)T  ~  h^^T-1).  We  then  find  the  scaling  p  ~  s-1/(y-1)  ~ 
Expressing  these  results  in  terms  of  an  equivalent  footpoint 
displacement  5  s  and  an  equivalent  B^  s  (2p)^'^^  as  suggested  by 

eqs.  (14a)  and  (14b),  we  find  B^  ~  Again,  these  scalings  are  in 

qualitative  agreement  with  the  results  of  Fig.  12a. 

The  results  presented  in  Figs.  5-11  and  Fig.  12b  have  a  bearing  on  the 
development  of  a  so-called  "open  field  configuration".  This  open  field 
configuration  is  defined  as  the  potential  solution  having  the  same  boundary 
condition  on  the  flux  <|<(x,o)  but  necessarily  having  a  sheet  current  along 
=  It  has  been  stated  (Barnes  and  Sturrock  1972;  Yang,  Sturrock,  and 

Antiochos  1986;  Aly  1984,1985)  that  the  field  approaches  such  a 
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configuration  as  the  footpoint  displacement  increases.  The  current  is 
conjectured  to  become  entirely  concentrated  into  the  sheet.  As  noted  in 
the  discussion  of  Fig.  7c,  the  transverse  components  of  the  current  density 
do  not  become  peaked  in  space  as  increases.  Also  the  magnitude  of  the 

transverse  current  density  jp  scales  as  Vax  ''ith  increasing  d,„gjj  and  thus 
increases  slowly  for  large  Furthermore,  the  total  transverse  current 

is  proportional  to  found  to  decrease  as  d^g^j^  is  increased. 

The  other  component  of  current  density  jz  =  becomes  peaked  near  \l> 

=  g/Q  (see  the  discussion  of  Fig.  7c).  However,  the  peak  of  jg,  at  x  =  a/2 
increases  slowly  for  large  dp,g^^  (Fig.  10).  (The  slow  increase  in  the  peak 
value  of  jg  for  d^gj^  ^  0.5  is  again  attributable  to  the  influence  of  the 

walls.)  Furthermore,  the  total  jg  current  in  the  peak  region  decreases. 
On  the  other  hand,  the  total  z-current  increases  at  a  considerably 

faster  rate  than  the  peak  current  density.  We  conclude  that  the  peak  in 
near  v|/  =  \|/q  contributes  a  rapidly  decreasing  fraction  of  the  total  current 
Ig  as  djjigjj  increases,  and  therefore  sheet  currents  with  progressively 
narrower  current  profiles  carrying  the  bulk  of  the  current  do  not  occur. 
The  increase  in  appears  to  be  due  to  the  same  cause  as  the  increase  of 
Ig  for  eq.  (8)  for  which  «  X^/(n^/a^  -  namely  that  the 

equilibrium  continues  to  expand,  with  nearly  uniform  current  density 
filling  the  space.  We  conclude  that  increasing  the  footpoint  displacement 
does  indeed  cause  the  flux  surfaces  to  expand  but  that  the  resulting 
equilibrium  approaches  a  configuration  with  finite  current  density  (j^)  and 
increasing  total  current  but  no  concentration  of  the  current  into 

sheets . 

In  the  present  paper,  we  have  adopted  the  conventional  approach  that 
the  footpoints  of  magnetic  field  lines  are  tied  to  the  infinitely 
conducting  and  massive  photosphere  and  have  found  no  bifurcations  (i.e.,  no 
loss  of  equilibrium)  with  respect  to  physically  specifiable  quantities.  In 
this  scenario,  all  the  flux  associated  with  the  arcade  is  contained  in  the 
corona.  If  the  photosphere  is  assumed  to  have  finite  conductivity  and 
finite  mass,  then  the  fluxes  can  extend  below.  In  a  recent  paper,  Chen 
(1989)  studied  the  behavior  of  "toroidal"  current  loops  in  which  the 
current  was  assumed  to  be  closed  in  or  below  the  photosphere.  In  this 
work,  a  circuit  parameter  e  s  ip/if  =  hp/L-p  was  found  to  play  a  role  in 
determining  the  stability  behavior  of  a  loop,  where  fp  (Lp)  is  the  flux 
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(inductance)  in  the  corona  enclosed  by  the  loop  and  (L>j)  is  the  total 
flux  (inductance)  including  the  flux  structure  below  the  photosphere 
associated  with  the  loop.  Thus,  the  quantity  e  parametrizes  the 
subphotospheric  flux/current  structure  relative  to  what  is  above  the 
photosphere.  In  this  model,  the  infinitely  conducting  and  massive 
photosphere  assumption  corresponds  to  the  e  =  1  limit.  In  this  limit, 
loops  are  found  to  be  stable.  If,  on  the  other  hand,  the  subphotospheric 
flux  is  large  enough  (e  less  than  some  critical  value  s^r^’  loop  can 
be  unstable  to  major  radial  expansion,  leading  to  a  wide  range  of  motion 
and  magnetic  energy  dissipation.  In  the  present  linear  arcade  model,  too, 
it  is  possible  that  if  the  fluxes  are  allowed  to  extend  below  the 
photosphere,  loss  of  equilibrium  or  instabilities  may  occur. 
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V.  Summary 


A  study  of  two-dimensional  MHD  equilibrium  of  linear  coronal  arcades 
has  been  presented.  An  iterative  method  has  been  used  to  solve  the  Grad- 
Shafranov  equation  with  prescribed  footpoint  displacement  d(il/)  or  entropy 
s(ij»).  This  method  involves,  as  an  intermediate  mathematical  step,  the 
prescription  of  and  p('J/).  A  modification  of  the  method  has  been 
developed  to  compute  equilibria  through  purely  mathematical  bifurcation 
points  which  can  occur  when  axial  field  and  pressure  are  prescribed  in 
the  intermediate  steps.  Sequences  of  equilibria  have  been  computed  with 
d(ij;)  or  s(^)  specified.  No  multiple  solutions  exist  if  d(\(<)  is  specified, 
consistent  with  the  suggestion  of  Jockers  (1978)  and  in  agreement  with  the 
recent  result  of  Zwingmann  (1987).  We  have  also  shown  that  no  multiple 
solutions  exist  if  the  entropy  s(i<()  is  specified.  The  entropy,  like  the 
pressure,  is  a  measure  of  the  internal  energy  of  the  plasma.  Generally,  if 
the  entropy  is  specified,  then  the  pressure  cannot  be  specified  and  vice 
versa.  In  the  limit  in  which  the  arcade  plasma  can  be  assumed  to  be 
adiabatic,  then  the  entropy  rather  than  the  pressure  is  physically 
specifiable.  In  this  case,  the  pressure-based  bifurcations  do  not  indicate 
the  possibility  of  an  eruptive  process.  We  have  described  physical 
conditions  under  which  the  entropy  may  be  specified,  and  discussed  the 
difficulty  in  posing  a  physically  meaningful  situation  in  which  the 
pressure  is  specifiable.  Our  results  indicate  that  increasing  footpoint 
displacement  and  entropy  increases  the  magnetic  energy  of  an  arcade  but  not 
the  magnetic  free-energy  to  drive  eruptive  processes  in  the  corona. 

We  have  also  investigated  the  formation  of  the  so-called  open  field 
configuration  and  formation  of  sheet  current  (e.g.,  Barnes  and  Sturrock 
1972;  Aly  1984,1985;  Yang,  Sturrock  and  Antiochos  1986)  and  found  that  the 
profile  of  the  current  density  indeed  becomes  more  peaked  with 
increasing  footpoint  displacement  but  that  the  total  current  in  the  peak 
region  decreases  as  the  total  current  I^  increases,  so  that  a  rapidly 
decreasing  fraction  of  current  is  carried  by  the  current  peak  at  the  center 
(x  =  a/2).  Based  on  our  solutions,  especially  Figs.  10,11  and  12(b),  we 
expect  no  sheet  current  formation  (i.e.,  concentration  of  current  at  x  ^ 
a/2)  to  takes  place  as  a  result  of  increasing  footpoint  displacement. 
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Appendix  A 


As  an  aside,  we  note  that  Clebsch  variables  have  been  used  by  other 
authors  in  specifying  footpoint  displacements.  That  is,  the  total  field  is 
written  as 


B  =  7a  X  73, 

where  a  =  \J;  and  3  =  z  -  f(x,y).  Writing  f(x,y)  =  0(x,y)dQ(«J;)  and  using  eq. 
(la),  we  find 

B-70  =  B^(M>)/d^(>p) 


Pds 

°  d^(^)  J  B 


with  an  indefinite  integral  along  the  field  line.  If  0  is  required  to 
increase  by  unity  in  integrating  from  one  footpoint  to  the  other,  we 
obtain,  from  eq.  (6a),  that  dQ(iJ^)  =  d(i)»).  Therefore,  specification  of 
f(x,0)  at  one  footpoint  (e.g.  where  >  0)  is  equivalent  to  specifying 
d('l') . 
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Appendix  B 


The  computational  method  we  employ  is  to  start  with  a  trial  For 

brevity,  we  begin  with  the  case  p  =  0.  We  iterate  eq.  (5a)  by  solving 


where  RC'J/)  =  -  B2('J')dB2/d<J;.  Given 
solve  eq.  (Bl)  for  'I'k+j.  Formally, 
It  is  helpful  to  introduce 
according  to  =  rL~l[R(\<^)]  +  (1 

the  inner  loop  iteration  'lv(s) 
properties  of  eq.  (5a)  with  the  trial 
out  the  iteration  past  the  mathemat: 
specifying  B^  as  an  intermediate  step, 
step  of  renormalizing  such  that 

fixed  during  successive  iterations. 

For  a  given  footpoint  displacement 
lines  to  compute  V'(\p),  defined  in  eq 
eq.  (6a)  with  a  relaxation  parameter  p 


=  R(^),  (Bl) 

a  standard  Poisson  solver  is  used  to 
this  can  be  expressed  as  =  L" 

a  relaxation  parameter  r  and  iterate 
r)<</j^^.  The  convergence  properties  of 
'<Tt+l(x)  depend  on  the  bifurcation 
function  B^f'J/).  In  order  to  carry 
cal  bifurcation  points  arising  from 
it  is  necessary  to  take  the  further 
the  total  axial  current  remains 

d(\j/),  we  integrate  along  the  field 
(6b).  We  then  compute  B2('J')  from 
such  that 


Bnew(^)  ^  pd(,|,)/V'(g»)  +  (1  -  p)B20^'^(g/).  (B2) 


The  method  of  inner-outer  loop  iteration  is  best  summarized  in  a  flow  chart 
as  follow.  First,  we  specify  the  footpoint  displacement  d(\|;).  We  then 

1)  guess  a  solution  to  eq.  (Bl)  by  giving  and  ^(x,y),  and  compute  the 

total  axial  current  I^; 

2)  iterate  eq.  (Bl)  with  a  (inner-loop)  relaxation  parameter  r.  If  is 

on  the  right  side  of  the  mathematical  bifurcations  (e.g.,  for  > 

0.35  for  the  B-series  equilibria  of  Fig.  8),  renormalize  B^f'J/)  at  each 
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iteration  so  that  I2  remains  unchanged  for  the  series  of  "inner  loop" 
iterations.  There  is  no  bifurcation  with  respect  to 

3)  compute  (eq.  I6bJ)  for  these  solutions  to  eq.  (5); 

4)  update  by  eq.  (B2)  with  a  second  (outer-loop)  relaxation  parameter 

p,  and  return  to  step  2. 

For  small  footpoint  displacement,  is  insensitive  to  changes  in 

footpoint  displacement  because  is  the  potential  solution  to  lowest  order. 
In  this  case,  the  outer  loop  iteration  (B2)  (with  relaxation  parameter  p  = 
1)  converges  rapidly.  For  larger  footpoint  displacement,  becomes 

quite  sensitive,  but  in  practice  it  is  found  that  the  iteration  oscillates 
if  it  diverges,  and  therefore  a  positive  relaxation  parameter  p  (with  p 
decreasing  with  guarantees  convergence,  provided  the  iteration 

process  is  properly  chosen  if  mathematical  bifurcation  with  B^  specified 
occurs  as  described  above.  This  is  discussed  further  in  Sec.  lie  for 
specific  examples.  This  method  of  inner  and  outer  loop  iterations 
generally  converges  very  rapidly  if  the  correct  relaxation  parameters  are 
chosen.  However,  it  often  requires  an  accurate  first  guess  for  <|»(x,y)  if 
footpoint  displacement  is  large.  Therefore,  we  generally  find  solutions  by 
varying  parameters  from  one  run  to  the  next,  starting  with  small  footpoint 
displacement . 

A  similar  scheme  is  used  to  include  pressure  when  the  entropy  s(\|/)  is 
prescribed.  The  analogue  to  eq.  (Bl)  is  R('J')  =  -dp/d\J/  and  p(<|/)  is 
renormalized  to  keep  I^  fixed  during  inner  loop  iterations,  as  discussed  in 
Sec.  IV.  Using  eq.  (6c),  the  analog  to  eq.  (B2)  is 

pnew(^)  ^  prs(^)/V'(<p)r  +  (1  -  p' )p^'l'^(«J/) .  (B3) 

Since  the  right  hand  side  of  the  G-S  equation  (5a)  is  -(d/di);)(p(\/;)  + 
B2('<')^/2],  the  prescriptions  for  computing  B2('J/)  and  p('l'),  and  their 
effects  on  the  final  solutions,  are  quite  similar.  This  is  most  easily 
seen  if  we  rewrite  eqs.  (6a)  and  (6c)  as 
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and 


i  =  I  d2(i|.)/V'(0»)^  (B4) 

pC'i/)  =  s(^0)/V'(^)y.  (B5) 


In  fact,  if  it  were  not  for  the  rather  insignificant  difference  between  y  = 
5/3  and  2,  there  would  be  a  formal  identification  (1/2)B2('|>)^  <-->  PC'!*)  and 
(l/2)d(\J;)^  s(t|/).  [Similar  statements  hold  for  the  isothermal  plasma,  in 
which  the  analogue  of  eq.  (B5)  is  y  =  1  with  s(\(>)  replaced  by  M'(\p).J 
Because  of  this  similarity,  the  iteration  scheme  eq.  (B3)  has  essentially 
identical  properties  as  eq.  (B2),  but  with  somewhat  different  optimum 
relaxation  parameters  in  the  outer  loop. 
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z-direction  and  a  '‘transverse”  peipendicular  to  z. 


FOOTPOINT  DISPLACEMENT  d 


Fig.  2  (Continued)  —  A  force-free  arcade  for  d(^)  given  by  Eq.  (16)  and  ^(x,0)  =  sin(irj:/a).  A 
relatively  small  amount  of  potential  flux  =  0.195  surrounds  the  current-carrying  plasma,  (a)  Flux 
(^)  surfaces,  (b)  Footpoint  displacement  versus  x.  (c)  Bj  versus  (d)  X  =  dBi/d  ^ 
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Fig.  2  (Continued)  —  A  force-free  arcade  for  given  by  Eq.  (16)  and  ypixfi)  =  sin(xj;/a).  A 
relatively  small  amount  of  potential  flux  \l/y  =  0.195  surrounds  the  current-carrying  plasma,  (a)  Flux 
(^)  surfaces,  (b)  Footpoint  displacement  versus  x .  (c)  versus  \f>.  (d)  \  =  dB^/d  ^ 


1.0  0.5  0 


'I' 

Fig.  2  (Continued)  —  A  force-free  arcade  for  d(^)  given  by  Eq.  (16)  and  ^(jc.O)  =  sin(TJc/a).  A 
relatively  small  amount  of  potential  flux  V'v  =  0-195  surrounds  the  current-carrying  plasma,  (a)  Flux 
(^)  surfaces,  (b)  Footpoint  displacement  versus  jt .  (c)  versus  (d)  X  =  dB^/d  ^ 
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Fig.  3  —  versus  (U...  Curve  1  corresponds  to  equilibrium  given  by  Eq.  (8)  and  curve  2  to 

equilibrium  given  by  Eqs.  (17)  and  (18).  Curve  3  describes  the  numerical  equilibrium  with  parame¬ 
ters  given  as  in  Fig.  2. 


max 

Fig.  4  —  (BJniax  and  the  total  current  as  functions  of  for  the  equilibrium  analytically  given  in 
Eqs.  (19)  and  (20)  with  v  =  1. 


Fig.  5  —  Numerical  equilibrium  with  boundary  condition  given  by  Eq.  (19)  for  k  =  7  and  with  foot- 
point  displacement  given  by  Eq.  (21).  The  current-carrying  region  is  surrounded  by  a  relatively  large 
amount  of  flux  =  0.25  (jCy  =  0.31),  W  =  0.40,  and  =  0.36.  (a)  Flux  surfaces,  (b)  Fooqx>int 
displacement  versus  x.  (c)  versus  tp.  (d)  X  versus  \p. 


Fig.  5  (Continued)  —  Numerical  equilibrium  with  boundary  condition  given  by  Eq.  (19)  for  =  7 
and  with  footpoint  displacement  given  by  Eq.  (21).  The  current-carrying  region  is  surrounded  by  a 
relatively  large  amount  of  flux  V-v  =  0-25  (Xy  =  0.31),  W  =  0.40,  and  d^a*  =  0.36.  (a)  Flux  surfaces, 
(b)  Footpoint  displacement  versus  x.  (c)  Bj  versus  \l/.  (d)  X  versus 


Fig.  5  (Continued)  —  Numerical  equilibrium  with  boundary  condition  given  by  Eq.  (19)  for  »>  =  7 
and  with  footpoint  displacement  given  by  Eq.  (21).  The  current-carrying  region  is  surrounded  by  a 
relatively  large  amount  of  flux  =  0.25  (jCy  =  0.31),  W  =  0.40,  and  (!„„  =  0.36.  (a)  Flux  surfaces, 
(b)  Footpoint  displacement  versus  x.  (c)  versus  (d)  X  versus 


1.0 
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0.5 

Fig.  5  (Continued)  —  Numerical  equilibrium  with  boundary  condition  given  by  Eq.  (19)  for  «'  =  7 
and  with  footpoint  displacement  given  by  Eq.  (21).  The  current-carrying  region  is  surrounded  by  a 
relatively  large  amount  of  flux  =  0.25  {Xy  =  0.31),  W  =  0.40,  and  =  0.36.  (a)  Flux  surfaces, 
(b)  Footpoint  displacement  versus  x.  (c)  Bj  versus  (d)  X  versus 
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1.0 


0.5 
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Fig.  6  (Continued)  —  Equilibrium  with  parameters  as  in  Fig.  5  except  dm„  =  0.61. 
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^p/  gp  = 


1.0  0.5  0 

Fig.  7  (Continued)  —  Equilibrium  with  parameters  as  in  Fig.  5  except  dm**  =  1-35. 
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Fig.  8. 
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